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ABSTRACT 

Extensive Monte Carlo modeling of a coded aperture 
y-ray telescope based on a high resolution liquid xenon 
TPC has been performed. Results on efficiency, back- 
ground reduction capability and source flux sensitivity are 
presented. We discuss in particular the development of a 
reconstruction algorithm for events with multiple interac- 
tion points. From the energy and spatial information, the 
kinematics cf Compton scattering is used to identify and 
reduce background events, as well as to improve the detec- 
tor response in the few MeV region. Assuming a spatial 
resolution of 1 mm RMS and an energy resolution of 4.5% 
FWIIM at 1 MeV, the algorithm is capable of reducing 
by an order of magnitude the background rate expected at 
balloon altitude, thus significantly improving the telescope 
sensitivity. 

1. INTRODUCTION 

Significant progress in experimental y-ray astrophysics 
can only be achieved with the development of instruments 
with improved imaging capability. Efforts to improve the 
sensitivity of a y-ray telescope without a corresponding im- 
provement in angular resolution will result in source confu- 
sion in the field-of-view (FO V) and will limit the mapping 
of extended source distributions. Two important scien- 
tific objectives which typify the need for detectors with 
improved imaging capability, along with good sensitivity 
and energy resolution, are the 1.8 MeV line emission from 
radioactive 26 A1 and the 0.511 MeV positron-electron an- 
nihilation emission. For both cases a precise study of the 
spatial distribution of the emission will uniquely identify 
the sources of radiation. Among the proposed novel tech- 
niques for imaging astrophysical y-ray sources in the low 
to medium energy region, the liquid xenon Time Projec- 
tion Chamber (LXe-TPC) is recognized as very promising 
and worth a strong R&D effort. Like an electronic bub- 
ble chamber, a TPC is capable of visualizing the complex 
histories of y-ray events with multiple interactions, ini- 
tiated by either Compton scattering or pair production. 
As a result of this imaging, efficient background rejection 
is also achieved, reducing the requirement for a massive 
anti-coincidence shield of the type that is required for ger- 
manium and sodium iodide detectors. From the energy 
and spatial information available for each point of inter- 
action, Compton kinematics enables the reconstruction of 


the direction of the incoming y-ray The direction is not 
unique unless the Compton scattered electron is also im- 
aged. Given the short range of low energy electrons in 
liquid xenon, this is not possible for a practical spatial 
resolution of the order of millimeters. In the LXe-TPC 
Compton telescope, as proposed by us in [1], the angu- 
lar resolution for the MeV region is limited by the small 
separation between two successive y-ray interactions. 

To overcome this limitation and release the strong re- 
quirement on the TPC spatial resolution, we propose to 
combine the imaging capability of the LXe-TPC with that 
of the coded aperture technique, for studying astrophysi- 
cal sources in the few MeV region with areminute accuracy. 
In a coded aperture telescope, a mask consisting of a two- 
dimensional array of opaque and transparent elements is 
placed between the source and a position sensitive y-ray 
detector. Every source within the FOV casts a shadow of 
part of the mask onto the detector. By properly decoding 
the pattern obtained at the detector plane, an image of 
the source is inferred. The successful application of coded 
aperture imaging in the low energy y-ray region has been 
demonstrated by two balloon borne telescopes : the New 
Hampshire "Directional Gamma-Ray Telescope”, DGT [2] 
and the Caltech “Gamma-Ray Imaging Payload”, GRIP 
[3], and more recently by the SIGMA telescope aboard the 
GRAN AT satellite [4]. All these telescopes use scintilla- 
tors as position sensitive y-ray detectors. The advantage 
of the LXe-TPC/coded telescope is in its unique capabil- 
ity to suppress background over a wide energy range, thus 
improving the sensitivity to weak sources, while maintain- 
ing an excellent detection efficiency, energy and angular 
resolution. 

In this paper the technical aspects of the LXe-TPC 
development are not considered. For this we refer to pre- 
vious publications by Aprile et al. [5]. We present here 
initial results from a Monte Carlo analysis of the expected 
performance of the LXe-TPC coupled to a coded aperture 
as y-ray telescope for 0.3-10 MeV. 

2. MONTE CARLO ANALYSIS AND RESULTS 

2.1 Telescope Model 

The model of the y-ray telescope used in the Monte 
Carlo simulation is schematically shown in Fig.l. It con- 
sists of a coded aperture mask, located 1 meter above the 
LXe-TPC y-ray detector. An active Csl shield has been 
included. In practice, the background rejection capability 



of the TPC reduces the need for massive anti-coincidence 
shield. The sensitive area of the TPC is 39 x 28 cm 2 . The 
total depth of liquid xenon is 10 cm. The detector’s dimen- 
sions and mask-detector separation have been determined 
on the basis of the existing gondola for the DGT telescope 
that we plan to use for the first balloon flight of the tele- 
scope [6]. The TPC is operated in the ionization mode. 
Both the energy and the spatial distribution of each ion- 
izing event within its sensitive volume is measured. The 
coordinates in the X-Y plane are inferred from the signals 
induced on two orthogonal wire planes, while the coordi- 
nate along the Z-direction is inferred from the drift time 
measured with respect to a zero time provided by the fast 
scintillation signal. As indicated by our first experimental 
results on the induction signals in a 3.5 liter LXe-TPC, to 
be published elsewhere, the spatial resolution in the plane 
is better than S/Vl2 where S is the spacing between the 
induction wires. For the spatial resolution on the coor- 
dinate along the drift axis we have previously measured 
o't — 0.18 mm [5]. For the Monte Carlo calculations, a 
resolution of <r r = a y = 1 mm and <r z = 0.2 mm has been 
assumed, along with the experimental energy resolution of 
4.5% FWHM at 1 MeV, extrapolated according to a E~ 1 / 2 
law [7]. 



Fig. 1. Schematic of the LX^TPC/coded mask 7 -m y telescope 


The mask is a 2 x 2 mosaic of a basic Uniformly Re- 
dundant Array (URA) [8]. It consists of a 43 x 41 ele- 
ment pattern of 0.91 cm x 0.68 cm x 1.2 cm thick tungsten 
blocks. With the assumed separation of 1 meter between 
mask and the detector, an angular pixel size of 0.65° x 
0.49* and a FOV of28 a x 20° is defined. Much finer angu- 
lar resolution is obtained by increasing the mask-detector 
separation. The precision in point source localization de- 
pends, however, also on the detector’s spatial resolution 
and the statistical significance of the source data. We have 
estimated that a precision at the arc minute level is possible 
with a practical spatial resolution of 1 or 2 nun. 


2.2 Reconstruction Algorithm based on Compton Kine- 
matics 

The interaction of 7-rays in liquid xenon is dominated 
by Compton scattering and pair production, for energies 
larger than few hundred keV. The typical event registered 
by the TPC is therefore characterized by multiple inter- 
action points before the original 7-ray energy is fully de- 
posited in the liquid. This is illustrated in Fig. 2 where the 
distribution of the number of interaction points for three 
different 7-ray energies is plotted. The probability that the 
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first interaction point is a Compton scattering is 90% at 
1 MeV, decreasing to about 25% at 10 MeV. To improve 
the signal/noise level in the few MeV region where multiple 
Compton events dominate, an algorithm based on the kine- 
matics of the Compton process was developed and tested. 
The algorithm uses the capability of the LXe-TPC to vi- 
sualize each incoming 7-ray event and to reconstruct the 
initial direction from the energy deposited in each inter- 
action point and the coordinates of the points. Consider 
an event in which a 7-ray produces N energy deposition 
points. Let Wo be the energy of the incoming photon, W, 
the energy of the scattered photon, Ei the energy of the 
scattered electron and cos 0,- the cosine of the scattering 
angle. If the coordinates and the energy depositions for 
each interaction poiut, as well as the right order of the 
successive interactions are known, then from energy and 
momentum conservation we can write for each interaction 
point: 


1 — cos 0,- = 


Wf 


1 

W.-i 


( 1 ) 


where W, _ W'i+j + Ei+ j, for i = l,Af - 1. For any 
given distribution of interaction points generated inside the 
sensitive volume, we now assume that the last interaction 
point is a photo-absorption and that the initial direction 
of the incoming 7-ray is known. Calling 


Ei 

1 — cos 0,- a " 

(1) becomes 

Wi -1 ■ Wi = a,. 

Similarly, for the (»' + 1)‘ A point, 


( 2 ) 

( 3 ) 


Wi ■ W i+1 = o, +l 


( 4 ) 


From (3) and (4), 


Wi = 


<*i ~ <*»+! 

Ei + Ei+i 


Also, if the total energy is contained, W 0 is 


( 5 ) 


2.3 Detection Efficiency 

The Monte Carlo calculations were based on the EGS4 
computer code [9]. The excellent gamma-ray detection ef- 
ficiency of liquid xenon over a wide energy range is shown 
in Fig.3. The full energy peak efficiency was calculated 
for gamma-rays in the range 0.3-10 MeV, normally inci- 
dent on the forward face of a 10 cm deep liquid xenon 
detector. The energy of the events in the sample was re- 
quired to be totally contained. The use of a coded mask 



reduces the detection efficiency by a factor of two, assum- 
ing 100% opacity of the closed mask elements. This is 
shown in the upper curve of Fig. 4. The effect of 3 mm 
stainless steel detector wall was also accounted for. The 
effect of the Compton reconstruction algorithm, described 
in section 2.2, on the telescope efficiency is also shown in 
Fig. 4 with the points □. Because the algorithm is tai- 
lored for events which start with a Compton scattering, all 
events with a single photo-absorption or which start with 
a pair production are rejected. To recover the detection 


W 0 = 


«i 

E? + E 3 + • • • Efj 


( 6 ) 


Alternatively, from energy conservation alone we have an 
independent estimate of the energy of the scattered photon 
at each point. We express this as W-. 


W{ — Ei+i + • • • E n (7) 

f° r * — 0,7V— 1. To test the validity of the assumption 
that the total energy is contained, as well as if the assumed 
initial direction of the 7-ray is kinematically possible, the 
function 

N-l 

/=Ew-^) 5 - (8) 

•=o 

is nunimized for all possible combination of points. The 
calculation takes into account the errors on the energy and 
spatial coordinates. 




efficiency at low and high energy, different reconstruction 
techniques can be applied. For example, the majority of 
low energy photons coming from the source will interact in 
the top region of the LXe-TPC. Fiducial volume cuts for 
single site events can therefore drastically improve the sig- 
nal/noise ratio. This is shown in Fig. 4 with the points •. 
All single-site events reconstructed outside the first 1 cm 
layer of liquid xenon from the top were rejected. For high 
energy events which start with a pair production, tracking 
of the elec tron/ positron can be used, on an event-by-event 
basis, to estimate the 7 -ray direction from the opening an- 
gle of the pair. Such an algorithm has however not been 
applied as one can see from the drastic drop in efficiency 
at energies above ~ 2 MeV. 

The capability of the reconstruction algorithm to re- 
ject events under the Compton continuum is shown in Fig. 
5, which shows the ratio of the total counts under the peak 
to that under the Compton continuum. A large improve- 
ment in this ratio is achieved for the MeV region. The 
excellent capability of the same algorithm to recognize 7 - 
ray events entering the TPC normally from the bottom 
was also tested. Results are shown in Fig. 6 . For this 
calculation total energy containment was not required. 



Fig 5. Ratio of events under the peak to that under the ComDton 
ccotmuum for the LXe-TPC/coded mask before and a£?vS 
reconstruction. 


2.4 Background Reduction and Flux Sensitivity. 

The ultimate flux sensitivity is perhaps the most rel- 
evant “figure-of-merit” of a 7 -ray telescope. Among the 
(actors which determine the sensitivity are: (i) the energy 
of the 7 -rays and (ii) the detector’s effective area, energy 
and spatial resolution, the exposure time and the total 
background rate. To calculate the background expected in 
the LXe-TPC/coded mask telescope at balloon altitudes, 
we have taken into account the dominant atmospheric and 
cosmic diffuse components, entering the forward aperture 
of the telescope or leaking through a thick (5 cm) Csl ac- 
tive shield. The flux and angular distribution of the at- 
mospheric 7 -rays used in the calculation were taken from 
the parameterized forms given by Costa et aL {10] and 
the cosmic diffuse spectrum used was that given in refer- 


ence [l l]. The internal backgrounds from natural radioac- 
tivity, cosmic ray induced radioactivity and activation of 
instrument materials has been neglected, as we expect to 
reject the majority of these single site events. Pulse shape 



discrimination of the scintillation light pulses offers ad- 
ditional background rejection capability, as it is strongly 
dependent on the ionization density of the event. For ex- 
ample, neutron induced events can easily be discriminated 
against, based on the light signature. Figure 7 shows the 
calculated background flux, before and after applying the 
Compton event reconstruction of section 2.2 and the fidu- 
cial volume cut. We assumed a float altitude of 3 g/ cm 2 , 
over Palestine, Texas. The integrated flux over the 0.1-10 
MeV region gives about 340 counts/sec, consistent with 
typical background rates measured at these altitudes. A 
background reduction of almost a factor of 10 is obtained 
by identifying gamma-rays which kinematically couldn’t 
have come through the FOV of the telescope. It is clear 
that even better background suppression can be achieved, 








on an event-by-event basis, using different algorithms for 
different gamma-ray energies, as previously discussed. 

From the calculated efficiency and background rate 
we have estimated the 3 <r minimum flux sensitivity of the 
proposed telescope. The results are shown in Fig. 8a and 
8b. For a typical exposure of 3 x 10 4 s, we find a value 
of 9 x 10~ 5 photons/cm 2 sec for 511 keV and of 6 x 10” 5 


Line Flux Sensitivity (3cr) 



Continuum Flux Sensitivity (3cr) 



Fig. 8 . The 3<r minimum (a) line and (b) continuum flux •ensi- 
tivity of the LXe-TPC/coded mask T-ray telescope. 


photons/cm 2 sec for 1.8 MeV. A continuum sensitivity of 
3 x 10- 7 photons/cm 2 sec is obtained at 1 MeV. The sen- 
sitivity curves of the other instruments, shown for com- 
parison, have been taken from Winkler (1991) [12]. When 
combined with the arcminute source localization accuracy, 
the high sensitivity of the LXe-TPC telescope makes it 
competitive with many satellite instruments, even with the 
much shorter observation time available in a ballloon flight. 

4.CONCLUSION 

Monte Carlo modeling of a realistic LXe-TPC/coded 
mask telescope has demonstrated the unique capability of 
tills instrument to suppress background events as a direct 
consequence of the imaging properties of the LXe-TPC. By 
using a reconstruction algorithm for multiple site events, 
based on the kinematics of Compton scattering, we find 
that the expected background flux at balloon altitude is 
reduced by nearly an order of magnitude, resulting in su- 
perior line flux sensitivity in the MeV region, compared to 
that of existing scintillator-coded mask or proposed seg- 
mented germanium/coded mask 7-ray telescopes. It is 
clear that this is only a first step towards the application 
of the liquid xenon technique for precise imaging of cosmic 
7- ray sources. The full potential will be realized when a de- 
tector with large effective area, improved energy resolution 
and sub-millimeter spatial resolution becomes practically 
operational as a Compton/pair production 7-ray telescope. 

This work was supported by NASA grant NAGW- 
2013. 
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